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ABSTRACT 

We have upgraded IRCAL, the near- infrared science camera of the Lick Observatory adaptive optics 
system, to add a dual-channel imaging polarimetry mode. This mode uses an optically contacted 
YLF (LiYF4) WoUaston prism to provide simultaneous images in perpendicular linear polarizations, 
providing high resolution, high dynamic range polarimetry in the near infrared. We describe the 
design and construction of the polarimeter, discuss in detail the data reduction algorithms adopted, 
and evaluate the instrument's on-the-sky performance. The IRCAL polarimeter is capable of reducing 
the stellar PSF halo by about two orders of magnitude, thereby increasing contrast for studies of faint 
circumstellar dust-scattered light. We discuss the various factors that limit the achieved contrast, and 
present lessons applicable to future high contrast imaging polarimeters. 

Subject headings: instrumentation: adaptive optics, instrumentation: polarimeters, instrumenta- 
tion:infrared 



1. INTRODUCTION 

In recent years, there has been a surge of interest in 
combining polarimetry with high angular resolution ob- 
serving techniques, such as adaptive optics (AO). This 
convergence is motivated both by polarimetry's diagnos- 
tic potential for studies of circumstellar dust, and by its 
ability to overcome atmospheric speckle noise. 

For spatially resolved circumstellar disks, the polar- 
ization of dust-scattered light can be used to constrain 
the nature of the scattering bodies, provi ding insights 
into dust grain structure and processing (e.g.lSilber et al.l 
[2000t iLucas etall[200l iGraham et al.Tl2007h . The Hub- 
ble Space Telescope can provide high angular resolution 
polar imetry with very high contrast (jHines fc Schneideil 
l2006l ). but the larger diameters of ground-based tele- 
scopes yield superior angular resolution. As currently 
planned, the James Webb Space Telescope will not in- 
clude any polarimetry capabilities. Thus the highest an- 
gular resolution near-infrared polarimetry in the foresee- 
able future must be obtained with AO. 

A key challenge in AO observations is overcoming the 
presence of an extended, temporally- variable speckle halo 
in the point spread function (PSF), t hat limits sensitivity 
for faint material near b r ight stars (iRacine e t al. 1999, 
lAime fc Soummed l2004t iFitzgerald fc Graham ,2006i ). 
Because atmospheric turbulence does not polarize speck- 
les, while dust-scattered light is strongly polarized, differ- 
ential polarimetry can separate them to reach the funda- 
men tal photon noise li mit for detection of faint material 
(e.g. iKuhn et al.l[2Ciol . 
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Seeking to obtain high-contrast polarimetry of dust 
around young stars, wc developed an imaging polarime- 
try mode for IRCAL (Lloyd et al. 2000), the InfraRed 
CAmera of the Lick 3- m Shane telescope and its laser 
guide star AO system (iMax et al.l Il997t iBauman et al.l 
l2002f l. IRCAL's differential polarimetry mode uses a 
yittrium lithium fluoride (YLF) WoUaston prism ana- 
lyzer and rotating half-wave plate modulator to reduce 
speckle halos by two orders of magnitude compared to di- 
rect AO imaging with the same instrument. We have al- 
ready pubhshed several studies of Her big Ae/Be stars us- 
ing this instrument, which includ ed brief d escriptions of 
the in strument itself (jPerrin et a l. 2004a, bl: iPerrin et al.l 
[200l . This present paper describes in greater detail the 
design and performance of the IRCAL polarimeter, and 
the data reduction methods adopted. Our goals are both 
to provide a detailed instrumental reference for current 
observations, and to present lessons learned that may be 
relevant for future AO polarimeters. 

The structure of this paper is as follows: After briefly 
restating some polarimetry theory (f|2]), we describe the 
design and construction of the IRCAL polarimeter (^ 
and present the results of various engineering tests of its 
performance (Q. We then discuss our on-the-sky ob- 
servation and calibration methods (fj5]) followed by the 
data reduction pipehne (S}6]). In Sj7]we evaluate the in- 
strument's achieved contrast. We close with a few ex- 
amples of astronomical data taken with the polarimeter, 
followed by a brief discussion and conclusion (Jj8|). 

2. THEORY FOR HIGH CONTRAST POLARIMETRY 

2.1. A Review of Polarimetry Fundamentals 

iTinbergenI (|1996D and iKelleil (|2002D provide excel- 
lent introductions to astronomical polarimetry, while 
lAdamson et al.l ()2005r ) summarizes the recent state of the 
art. 
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We briefly repeat here a few basic results to estab- 
lish notation. The polarization of light is usually rep- 
resented by the Stoke s vector [J, Q,U,V] (IStokesI 11853 : 
IChandrasekhail Il946f ). The usual astronomical con- 
vention is for the +Q direction to be oriented north- 
south, and +U northeast-southwest, with angles increas- 
ing counterclockwise from north to east. Linear polariza- 
tion can also be expressed in terms of polarized intensity, 
P = \fCp + lP , and position angle 6* = 5 arctan([//(5). 
For astrophysical situations involving the scattering of 
light, circular polarization is usually (though not always) 
small compared to linear polarization, so we shall gener- 
ally drop Stokes V . The normalized polarized intensity 
P/I is referred to as the degree of polarization, polar- 
ization fraction or percent polarization. Notation is not 
always consistent: some authors use P to refer to polar- 
ized intensity while others use it for degree of polariza- 
tion. In this work, capital /, Q, ?7, V, P will always refer 
to intensities (e.g. with units of janskies or Jy arcsec"^), 
not normalized quantities. 

Visible and infrared astronomical detectors are rela- 
tively insensitive to polarization, so to measure polariza- 
tion we must encode it in variations of total intensity. 
The simplest method uses a rotatable linear polarizer 
to allow only one polarization to reach the detector. To 
fully recover the three unknowns (Stokes /, Q, U) requires 
measurements from at least three suitably-chosen angles. 
This method has the disadvantages that (a) half the light 
is thrown away, and (b) Stokes parameters are obtained 
from subtraction of non-simultaneous images, so atmo- 
spheric variations cause spurious apparent polarization. 
This is particularly a problem for AO observations, which 
have complex and time-variable PSFs. 

Dual channel polarimetry avoids this by obtaining 
simultaneous measurements of perpendicular polariza- 
tions. A WoUaston prism ( ii3.1|) placed in the coUimated 
beam of a conventional focal-reducer style camera pro- 
duces two perpendicularly-polarized images of a source. 
The resulting difference image should be immune to vari- 
ations in seeing or atmospheric transparency. A polar- 
ization modulator (such as a rotating half-wave plate) 
is required to measure both Stokes Q and U . Fur- 
ther modulations that swap polarization states between 
the two sides of the detector can be used to minimize 
the effects of flat flcldi ng and other i nstrumental errors 
(|Kuhn et a l. 2001; Pat at k, Romaniello 2006). Carefully 
designed dual-beam, non-AO polarimeters can be made 
remarkably robust against both instrumental and atmo- 
spheric systematics, resulting in sen sitivities to polar- 
ization frac tions as low as 10~^ (e.g. jKemp et allllQSTt 
iHough et a l. 2006). 

2.2. High Contrast Differential Imaging Polarimetry 

Before we proceed into the details of the IRCAL po- 
larimeter, we must flrst establish the instrument's sci- 
entiflc aims. Fundamentally, our goal is to measure 
spatially varying polarization from resolved circumstellar 
dust. No attempt is made to measure the total net po- 
larization of any unresolved source. An AO polarimeter 
is far from an optimal tool for measuring net polariza- 
tion, due to the complicated optical train and its inherent 
instrumental polarization. (While modulation can elim- 
inate the instrumental polarization effects of optics after 
the modulator, for AO systems it is typically impractical 



to locate the modulator before the entire AO optical sys- 
tem.) For point sources, much more accurate and precise 
measurements are possible using "classical" dual-beam, 
non-AO polarimeters incorporating on-axis optics with 
low instrumental polarization or modulator s located very 
early in the optical train (e.g. [Potter fc Graham ,20051 : 
IHough et al.ll2006l : lMasiero et alIl2007[K 

While such instruments can be made extremely robust 
against atmospheric and flat field effects, the situation is 
somewhat more complicated for AO imaging polarime- 
ters that attempt to resolve polarized structure hidden 
within the PSF halo. Detector limitations inevitably con- 
strain us to modulate more slowly than the atmospheric 
timescale. Any residual image motion not corrected by 
the AO system must be compensated for by registering 
images prior to subtraction. Because of this, imaging po- 
larimeters remain sensitive both to uncertainties in the 
fiat field and errors in the registration process. As we 
discuss below, these factors limit the achieved contrast 
of our instrument. 

Furthermore, while AO differential polarimetry does 
produce a polarized intensity image P that is robust 
against atmospheric errors, the total intensity image / 
is identical to that from regular AO imaging, including 
the PSF speckle halo. Thus P/I will be a strongly biased 
quantity, giving in practice only a lower limit to the true 
polarization fraction. This situation will improve with 
the next generation of AO polarime ters (e.g. those in 
the SPHERE and GP I instruments; iBeuzit et all 120061 : 
[Macintosh et al.l 12006) which will benefit from greatly 
improved AO PSF quality and coronagraphic suppres- 
sion of the PSF halo. 

3. DESIGN & IMPLEMENTATION 

3.1. General Design Considerations 

A WoUaston prism (|Wollastonlll82d) deflects perpen- 
dicular polarizations oppositely, with the angle of each 
output beam approximately equal to 

S ~ An tan a 

where a is the prism angle and An is the birefringence. 
If such a prism is placed conjugate to the telescope pupil, 
it will split the incident beam to create two polarized im- 
ages on the flnal detector plane. The Lagrange invariant 
implies that each of these images will have an apparent 
sky-projected deflection of 

Vojcctcd ^ I ) An tan a (1) 

V ^tcl ) 

where -Dpupii and Dtci are the pupil image and telescope 
diameters (5 mm and 3 m respectively in our case). The 
polarimetric field of view is maximized if the total sepa- 
ration 2(5pi.ojcctcd is half the detector's total field of view. 
For IRCAL, this field of view is 19.4", implying an op- 
timal (5projected ~ 4.9". Thus once we have chosen our 
material (and therefore An), equation [T] may be solved 
f or the n ecessa ry prism angle a. 

lOlival ()1997D proposed using a wedged double WoUas- 
ton for simultaneous measurement of all four polarization 
states, and thus simultaneous extraction of Stokes Q and 
\J . However, this approach requires dividing the pupil in 
half, resulting in a loss of half the angular resolving power 
of the system, so the wedged double WoUaston is not an 
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Fig. 1. — A simplified schematic of the IRCAL polarimeter, not to scale but indicating the order of the various optical components. The 
box represents IRCAL's vacuum dewar. Labeled in bold are the polarimetry components added to the original AO system and camera: the 
rotating waveplate, the aperture mask, and the WoUaston prism mounted in the first of IRCAL's two internal filter wheels. The WoUaston 
prism splits the coUimated beam, forming on the detector two simultaneous images with perpendicular polarizations. Rotation of the 
waveplate modulates the intensity of the two images, allowing the measurement of Stokes Q and U. 



optimal solution for high resolution imaging. The pupil 
division also means that the Q and U beams see different 
portions of the atmosphere and will have PSFs at least 
as divergent as if they were non-simultaneous. Further- 
more, obtaining a given signal-to-noise level in both Q 
and U requires just as much total exposure time as in 
a regular dual-beam polarimeter. We therefore chose to 
develop a classical single- WoUaston polarimeter. 

We also chose to use a conventional rotating half-wave 
plate as the polarization modulator, instead of alterna- 
tives such as liquid crystal variable retarders (LCVRs) or 
ferroelectric liquid crystals. Ferroelectric liquid crystals 
are unsuitable given their too-fast modulation compared 
to typical infrared detector readout times. LCVRs pro- 
vide spectacular levels of precision in solar polarimetry, 
but are generally not achromatic, are difhcult to obtain 
for wavelengths beyond 1.8 /xm, and have birefringence 
that is temperature-dependent. In contrast, half- wave 
plates are available with excellent achromaticity and low 
sensitivity to temperature variations. Since we desired a 
single modulator for the entire 1-2.5 /im band, we use an 
achromatic waveplate as our modulator. Figure [1] shows 
the overall layout of the IRCAL polarimeter including 
the above-mentioned components. 

3.2. Choice of material for the WoUaston prism 

iQliva et al.l ()1997t ) provides an overview of candidate 
birefringent materials for infrared WoUaston prisms. We 
present in Appendix A an updated list of such materials, 
along with recent references for their optical properties. 

For an imaging polarimeter, it is important to choose 
a material whose birefringence Art does not vary 
much within each operating bandpass. A wavelength- 
dependent An causes different beam deflections for each 
wavelength within a bandpass, elongating the final im- 
age. Minimizing this effect, called "lateral chromatism", 
is particularly important for a high contrast polarime- 
ter: not only does lateral chromatism decrease angular 
resolution, it can greatly degrade PSF subtraction qual- 
ity by oppositely distorting the two polarizations' PSFs. 
Lateral chromatism may be quantified over a given wave- 
length range Ai < A < A2 by Oliva's parameter 

AnAi - AnAs 

This parameter is the birefringent analogue to the Abbe 



number V, which measures dispersion in optical mate- 
rials. High values of Vau correspond to low amounts of 
birefringent dispersion. The image elongation over the 
range Ai — A2 will be e = S/Vau- By selecting a mate- 
rial to maximize Vau, we can minimize blurring due to 
lateral chromatism. 

We chose the unaxial birefringent material YLF (yt- 
trium lithium fluoride, LiYF4, pronounced "yilf" by its 
users in the laser industry) because its An is remark- 
ably constant with wavelength, causing negligible lateral 
chromatism. Typically, Vau > 100 is considered good, 
corresponding to chromatic elongation less than 0.05" for 
IRCAL, better than the diffraction limit. YLF exceeds 
this (spectacularly at H, with Vau — 1290) and is readily 
available due to its widespread use in the laser industry. 
Its refractive index, 1.45, is similar to that of standard 
optical components, simplifying the optical design and 
manufacturing, and its birefringence. An = 0.022, al- 
lows the necessary beam divergence to be obtained with 
an reasonable thickness of YLF, roughly 4 mm. 

Three similar WoUaston prisms were fabricated by 
Onyx Optics of Dublin, C A and were received at Berkeley 
in March 2002 (Figure [2]). The pieces were joined using 
optical contacting and both exterior surfaces received a 
1-2.5 /im antireflection coating. The 30° prism angle pro- 
duces a beam divergence of 1.45°, giving a sky-projected 
separation of 8.6" between the two beams (close to the 
optimal value of half the field of view) . All three prisms 
performed as expected in laboratory tests. Two of the 
prisms exhibited minor surface defects, so we selected the 
third for astronomical use. That WoUaston was installed 
in one of IRCAL's internal filter wheels, in coUimated 
space just before the cold pupil. 

A concern with any birefringent material is that the 
differing optic axes will have different thermal expansion 
coefficients. For WoUaston prisms constructed by bond- 
ing two surfaces with perpendicular axes, temperature 
changes create stress between the two faces. YLF has 
fairly similar expansion coefficients of 5 x 10^^ K^^ par- 
allel and 8 x 10~^ perpendicular to the fast axis. 
Our prism has tolerated repeated thermal cycling to liq- 
uid nitrogen temperatures with no ill effects apparent 
after five years. 

One caveat about YLF is it is relatively expensive and 
can be hard to obtain in large sizes, such as needed for 



Fig. 2. — Left: The IRCAL dewar mounted on the Lick AO optical bench at the Cassegrain focus of the 3-m Shane telescope, highlighting 
the waveplate mechanisms. A rotation stage allows modulation of the polarization using the waveplate, and a linear translation stage allows 
inserting/removing the waveplate from the beam. Right: Our YLF WoUaston prism, manufactured by Onyx Optics of Dublin, CA. This 
prism has been installed in one of IRCAL's two cryogenic filter wheels since March 2002. 



WoUastons on 8-10 m telescopes. The greater angular 
resolution of large telescopes places a stringent constraint 
on lateral chromatism, requiring high Va„ materials like 
YLF for diffraction-limited performance. However, a 
large YLF WoUast on was recently f abricated for the Sub- 
aru Telescope (Ho dapp et al.ll2"007f ). demonstrating that 
the crystal fabrication challenges can be overcome. 

3.3. Waveplate and Aperture Mask 

We obtained an achromatic waveplate from Karl Lam- 
brecht Corp., fabricated from MgF2 and crystal quartz. 
The waveplate is optimized for 1-1.6 /im, and provides 
0.500 ± 0.005 waves of retardance over that range. At 
longer wavelengths, the retardance differs from half a 
wave by about 0.05 waves averaged across the Kg band. 
As a result the polarizing efficiency of the instrument is 
degraded slightly at long wavelengths, as discussed be- 
low in §4. The waveplate is installed in a Newport PR50 
rotation stage, itself mounted on a translation stage for 
insertion into the optical path (Figure [5]). For practical 
reasons, the waveplate is located at the entrance win- 
dow to IRCAL. Ideally, it should have been placed as 
far upstream as possible, since modulation allows any 
instrumental polarization from post-modulator optics to 
be separated from the astronomical signal, but given the 
constraints of the AO bench that was not feasible. As 
a result, polarizing optics before the modulator, such as 
tilted mirrors and the dichroic beamsplitter, introduce a 
nonzero instrumental polarization. 

To prevent overlap of the two polarized images, we 
installed in IRCAL's cold aperture wheel a rectangular 
mask corresponding to half the field of view, essentially a 
9" wide slit. As fabricated, the mask is slightly oversized, 
resulting in an 0.4" (6 pixels or 240 /Ltm) beam overlap 
in the middle of the detector. This overlap region is 
discarded during data reduction. 

3.4. Limitations and Imperfections 

The IRCAL polarimeter's performance is limited by 
several factors in its design and construction. For in- 
stance, errors in flat fielding prevent completely accurate 
subtraction of the two polarizations. Modulation reduces 



the impact of this, but cannot eliminate the problem en- 
tirely. Even though both polarizations are alternately 
measured on both sides of the detector, the images must 
be shifted to register them before subtraction, which de- 
pends on accurate flat fielding to maximize the cancella- 
tion of starlight. 

A second factor that limits real-wo rld polarimeters i s 
registration and shifting errors (e.g. lApai et al.l[200l . 
IRCAL's coarse pixel sampling (which only Nyquist- 
samples the PSF at wavelengths > 2 /xm) exacerbates 
this. The undersampled nature of J- and _ff-band 
data limits the precision possible in registration and in- 
troduces aliasing that compromises image subtraction. 
Even were the sampling improved, accurately registering 
images would still be challenging, particularly for deep 
exposures saturated on the central core. 

The placement of the WoUaston prism in a filter wheel 
before the cold pupil results in the two polarized beams 
already having diverged by the time they reach that pupil 
stop. This "pupil shear" causes each channel to have 
a slightly different effective pupil and significantly in- 
creases the non-common-path aberrations between the 
two beams. Traced backward through the telescope and 
onto the sky, each channel sees a slightly different col- 
umn of partially-corrected atmosphere, displaced from 
the other by ~ 5% of the pupil diameter. This results in 
non-common-path atmospheric wavefront errors, which 
cannot be compensated for by waveplate modulation. It 
would have been much better to place the WoUaston af- 
ter the cold pupil, but the layout of the existing filter 
wheels didn't allow this. 

IRCAL's aperture mask wheel is not a precision mech- 
anism. Its positioning has low repeatability, and it oc- 
casionaly shifts unpredictably (apparently as a result of 
electrical crosstalk from one of the other wheel mecha- 
nisms). Shifts in the aperture mask change the field of 
view and flat fleld for the observations. The fllter wheel 
containing the WoUaston is less sensitive, but slight off- 
sets from night to night change the relative position of the 
two polarized images, which the data reduction pipeline 
must deal with. 

Compared to these factors, any imperfections in the 
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polarimetry optics themselves, and/or the finite preci- 
sion of the waveplate rotation (0.1°), are entirely neg- 
ligible. WoUaston prisms frequently achieve extinction 
ratios better than 10~^; at this time, we cannot more 
precisely quantify the performance of our prism but be- 
lieve that its extinction ratio and polarizing efficiency are 
excellent. 

4. COMMISSIONING TESTS 

The prism was installed into IRCAL during 2002 
March, and had first light on 2002 March 29 UT. Ob- 
servations then and on subsequent nights have allowed 
us to quantify the performance of the system. 

The polarimetry optics have negligible effect on the 
achieved AO image quality. They do contribute addi- 
tional non-common-path wavefront error between the sci- 
ence detector and wavefront sensor (including about 1.3 
mm of focus and 85 nm RMS of other aberrations, prin- 
cipally astigmatism) but by applying appropriate offsets 
to the deformable mirror, we can fully compensate for 
these aberrations. This results in internal Strehl ratios of 
0.92 ± 0.02, essentially identical to performance without 
the polarimetry optics^. The resulting orthogonal polar- 
ized PSFs appear very similar but not perfectly identi- 
cal, confirming the presence of small non-common-path 
wavefront errors between the two polarizations. 

We observed the Lick AO system's internal calibra- 
tion source to determine the throughput of the WoUaston 
prism and waveplate, including both reflective losses and 
absorption. At longer wavelengths, YLF becomes more 
transparent, but the waveplate's antireflection coating 
performs less well. This results in a net transmission 
through both optics of 92.3 ± 0.3% at J and H , decreas- 
ing to 88.5 ±0.5% at K^. 

To calibrate the polarizing efficiency, we placed in front 
of the AO calibration source a 1-2 /im IR Polaroid from 
Edmund Scientific, and observed its apparent brightness 
as a function of waveplate angle, for each of J, H, and 
Ks- We fit a model to these data using the Levenberg- 
Marquardt least-squares method to extract the polariz- 
ing efficiency. The polarizing efficiency measured in this 
way is 1.00 ± 0.005 for J and H, and 0.95 ± 0.01 for Ks 
band. This decreased polarizing efficiency at Kg is due 
entirely to the waveplate's retardance differing slightly 
from half a wave at wavelengths beyond 1.8 /^m. 

4.1. Instrumental Polarization 

Instrumental polarization of some degree is found in 
all optical systems, particularly those that depart from 
circular symmetry, such as by having oblique refiections. 
The Lick AO optical path involves many such reflections, 
and thus moderate instrumental polarization is expected. 
Accurately modeling instrumental polarization is chal- 
lenging, as oxide layers on mirrors an d other subt le im- 
perfections can play substantial roles ()Kelleij[2002f ). For 
very high contrast applications, the polarizing effects of 
thin metal films must be considered (jBorn fc WoUl[l959l : 
[Breckinridge fc: Oppenheimeij|2004) . 

^ Strehl ratios are measured using an IDL-based PSF fitting 
routine developed by Bruce Macintosh and incorporated into the 
IRCAL quicklook software. There are numerous pitfalls in accu- 
rately measur ing Strehl ratios, particularly for undcrsampled data 
IjRoberts et al.i 12004 ). so Strehl ratios quoted for IRCAL are un- 
certain by 5 — 10% in an absolute sense 



We observed unpolariz ed standard stars (|Gehrelslll974l : 
lAgeorges fc Walshlll999| ) to calibrate instrumental polar- 
ization. We measured the apparent polarization of each 
standard via aperture photometry on the reduced data, 
and interpreted the apparent polarizations as the amount 
of instrumental polarization. For all wavelengths, the de- 
rived instrumental polarizations are around 2.5% in the 
—Q direction. Due to the Cassegrain conflguration, these 
values are extremely stable with time, changing < 0.1% 
in observations of standard stars more than a year apart. 
Based on comparison of these measurements to a Mueller 
matrix model of the optical train, the instrumental po- 
larization appears to be primarily caused by the dichroic 
that separates the visible wavefront-sensing and IR sci- 
ence channels. This large contribution is unsurprising 
for a complex multilayer optic at oblique incidence. We 
did not attempt to measure circular instrumental polar- 
ization or crosstalk between circular and linear polariza- 
tions. 

The optimal method to remove instrumental polariza- 
tion depends on the science goals. We adopt an iterative 
optimizing subtraction process that simultaneously re- 
moves linear instrumental and interstellar polarization, 
described in tj6.1l 

4.2. Pixel Scale 

Lloyd et al. (2000) previously measured IRCAL's pixel 
scale as 0.0756" ± 0.0002 pixel"^. Subsequent obser- 
vations have allowed a refined measurement, revealing 
slightly different plate scales in R.A. and Dec. This indi- 
cates that IRCAL is anamorphic, a result not unexpected 
given that its internal off-axis parabolic mirrors are not 
of unit magnification. 

To quantify the anamorphism, we obtained LGS AO 
observations of the core of the globular cluster M 53. We 
extracted point source coordinates using the Star finder 
IDL package, and compared them to coordinates de- 
rived from HST W FPC2 imaging of the same cluster 
(jPiotto et al.ll200^ . We simultaneously fit the displace- 
ment, rotation, shear, and x- and y-magnification be- 
tween the two datasets using the MPFIT Levenberg- 
Marquardt nonlinear least squares optimizer^ . This anal- 
ysis yielded plate scales of 0.0754" pixel"^ in R.A. and 
0.0780" pixel-i in Dec, with statistical error ±0.0005" 
pixel^^. The detector's axes are orthogonal within our 
0.5° measurement uncertainty. 

To correct the anamorphism, we developed an IDL rou- 
tine ircal_dewarp, based on Keek's nirc2warp by An- 
tonin Bouchez. This routine transforms a distorted input 
image onto a uniform 0.040" grid via two-dimensional cu- 
bic spline interpolation. 

To assess the performance of this correction, and to 
measure any shifts in the absolute orientation of IRCAL 
over time, we routinely observe binary stars drawn from 
the Suggested Calibration Targets list of the W ashing- 
ton D ouble Stars (WDS) Orbits Catalog (Hartkopf et al.l 
120011 ). The observed binary separations are compared to 
values computed from the WDS orbital elements. Be- 
cause the WDS Orbits Catalog is itself imperfect (most 
of the binaries with suitable separations and bright- 
nesses have orbits graded 4 or 5, the lowest grades) some 
amount of discrepancy is expected. Since many entries 

* http://cow.physics.wisc.edu/ craigm/idl/fitting.html 
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in the WDS do not contain estimated uncertainties, we 
make no attempt here to separate out these two effects. 
A conservative estimate is that the pixel scale in post- 
anamorphism-correction IRCAL images is uncertain by 
no more than ^2.5% in linear size, and by ^1.0° in po- 
sition angle. 

The cubic spline interpolation employed is a good ap- 
proximation to the theoretically optimum sine interpola- 
tion, but not perfect, particularly given the non-Nyquist- 
sampled nature of short-wavelength IRCAL data. "Ring- 
ing" artifacts (the "Gibbs phenomenon") are some- 
times seen near stellar PSF peaks in dewarped im- 
ages. It is worth noting that the ircal_dewarp rou- 
tine is not strictly flux conserving, nor does it account 
for the small higher-order distortions present in IR- 
CAL (M. P. Fitzgerald, private communication). At- 
tempting to characterize and compensate for these ef- 
fects, either through heuristic algorithms such as Drizzle 
(|Fruchter fc HooM T20021 or through more formally cor- 
rect Fourier interlacing methods (Bracewell 1978, Lauer 
1999), is beyond the scope of this paper. 

5. OBSERVING TECHNIQUES AND ASTRONOMICAL 
CALIBRATION 

5.1. Observing Procedure 

Our standard observing procedure is as follows: Once 
the target has been acquired, the observer configures set- 
tings such as filter and exposure time. The camera's soft- 
ware automatically rotates the waveplate and integrates 
at each position. We use eight distinct waveplate posi- 
tions: 0, 22.5, 45, 67.5, 90, 112.5, 135, and 157°. This 
provides four redundant measurements for e ach of Stokes 
Q and U, reducing instrumental s ystematics ()Kuhn et al.l 
120011: iPatat fc Romaniellol [2001 . To obtain the neces- 
sary high dynamic range, for most targets it is neces- 
sary to combine both short and long exposures, typically 
between 0.5 and 90 s, depending on source magnitude. 
Small dithers after each set of eight exposures are used 
to reduce flat field effects and/or increase field of view. 

5.2. Flat fielding 

As flat-fielding errors are expected to be a major con- 
tributor to systematic limits on the detection of faint 
polarized circumstellar material, several methods were 
pursued for flat-fielding the polarimeter. The IRCAL 
aperture mask and filter wheel positioning mechanisms 
have frustratingly poor repeatability. Because of this, the 
aperture mask position can only be approximately repro- 
duced from night to night (typically within 5-10 pixels), 
so it is necessary to take flats again every night. Some 
authors have advocated computing flat field calibrations 
using observations of extended sources such as planets 
(|Kuhn et al.lll99ll ). but such targets are not always vis- 
ible. We instead opt for the more traditional sky as flat 
calibrator. 

Because the twilight sky is itself highly polarized, it is 
far from ideal for flat fielding a polarimeter. Dome flat 
screens are typically only slightly better. To get around 
this problem, polarized sky or dome exposures taken with 
different waveplate angles are averaged to make "pseudo- 
unpolarized" flat frames. This method only works per- 
fectly if the sky brightness and polarization are constant, 
a requiremen t notably not satisfied during twilight (see 
ICronin et al.l ()2006f ) for discussion of sky polarization 



during twilight and its variation with time). IRCAL's 
small pixel scale allows it to take unsaturated exposures 
on the daylit sky hours before sunset, at which point the 
sky polarization is more stable. 

From each night's set of "pseudo-unpolarized" frames 
we compute one master median flat per wavelength, 
which we use for all waveplate angles. Errors in the 
final master flats are estimated to be a few parts in a 
thousand. Some authors have reported slight benefits 
from using a flat fram e computed in dividually for each 
waveplate angle (.Ageorges fc Walshlll999 ). We did not 
investigate that approach in detail, but note that it may 
offer future potential for improvement. 

The nightly flat fields are also used to measure the 
location on the detector of the two polarized subimages, 
using a Hough-transform edge detection algorithm to find 
the edges of the illuminated fields of view. The detected 
masks are undersized by 5 pixels to avoid light scattered 
from the aperture mask's edges. 

5.3. Sky Subtraction 

At night the infrared sky background is primarily 
unpolarized, vibrationally excited emission from atmo- 
spheric molecules. In dual-channel differential polarime- 
try, this unpolarized background will subtract out of the 
Q and U images. Thus there would be no need to worry 
about sky subtraction, if our sole concern was the po- 
larized images. However, since we also want to mea- 
sure total intensity for each target, we perform the usual 
sky subtraction. We subtract a distinct sky frame for 
each waveplate angle, obtained either from separate sky 
frames or as a median of the dithered science frames with 
that angle. This process also allows for the subtraction 
of any polarized sky light such as scattered moonlight, 
but as mentioned above, this is usually negligible. 

5.4. Polarization Angle Calibration 

The Stokes reference frame must be oriented relative to 
some known reference angle. (Equivalently, we must de- 
termine the rotational zeropoint for the waveplate mech- 
anism.) We do this using the twilight sky, which is 
strongly linearly polarize d perpendicula r to the direction 
pointing toward the sun (jCronin et all [2006 ). Using the 
sky as a calibrator allows us to use all pixels of the de- 
tector and thus rapidly obtain high signal to noise. In 
practice, the same data are used for flat fielding the de- 
tector and for calibrating the position angle'^. 

5.5. Photometric Calibration 

Photometric stand ards (e.g. iPersson et al.l 119981 : 
iHawarden et al.|[200lD are observed using the same pro- 
cedure as science targets. After the standard polarime- 
try reduction (described in we measure their fluxes 
within a 7" diameter aperture, plus encircled energy 
curves for aperture correction. Adaptive optics photome- 
try is a notoriously problematic area due to PSF variabi l- 
ity fsee lFitzgerald fc Grahamll2006l : iSheehv et al.ll2006f) . 
Changes in sky conditions and guide star brightness in- 
evitably cause wavefront correction to vary between sci- 
ence target and calibration source, causing photometric 

^ We note in passing that [Harrington et aLI 112006!) report difficul- 
ties using the sky as a position angle calibrator in Hawaii, because 
of specular reflection of the sun by the ocean, but we have found 
that it works adequately for our purposes at Mt. Hamilton. 
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errors as the encircled energy curve changes. We expect 
an uncertainty of 5-10% in the resulting photometry in 
typical conditions. 

6. THE IRCAL POLARIMETRY DATA REDUCTION 
PIPELINE 

Data reduction for imaging polarimetr y has been dis- 
cussed b efore by a nuinber of authors (iGledhill et al. 
1991; Wo lstencroft et al.lll995HKuhn et al.l 120011: iPottCT 
2003al: |PataL& Romaniellol 120061). The reduction algo- 
rithms for IRCAL build on these earlier works but also 
feature a number of elements specific to IRCAL. We de- 
scribe here the data reduction process in detail, in the 
hope of sharing lessons from IRCAL with other aspiring 
differential polarimetrists. Furthermore, this is an appro- 
priate setting to air dirty laundry: that is, to highlight 
what doesn't work about IRCAL polarimetry as well as 
what does, and what aspects should be improved in fu- 
ture designs, such as the planned polarimetry mode of 
the Gemini Planet Imager. 

6.1. General Considerations 

Our fundamental goal is to use polarimetry to cancel 
out the unpolarized stellar PSF and reveal faint circum- 
stellar material. In practice, complications such as in- 
strumental and interstellar polarization mean that the 
star's PSF is not truly unpolarized. For distant or young 
and embedded sources, interstellar polarization can be 
substantial, at a level of a few percent. While it is 
possible to remove interstellar polarization based on ob- 
servations of nearby stars (Heiles 2000), this is fraught 
with uncertainties. Not least of these is the assumption 
that interstellar polarizations vary in a predictable and 
uniform way between stars- an assumption unlikely to 
be satisfied in the clumpy, dusty environments of young 
stars. Is there a better approach? 

We can instead allow the net polarization of the star 
(arising from both interstellar and instrumental polar- 
ization) to be a free parameter that we solve for. We 
do this by finding the scaling coefficients that minimize 
residual starligh t in | Q J and \U\, following the approach 
used by iPotten (|2003ah . This method refies on the as- 
sumption that the total flux is dominated by inherently 
unpolarized light, so that any net polarization seen in the 
data is the result of instrumental or interstellar polariza- 
tion. For our science targets, this is generally an good 
approximation, since scattered light is no more than a 
few percent of the total. 

This approach is robust even in the presence of scat- 
tered light. For a symmetric distribution of dust around 
a star, the polarization position angle varies around the 
disk and the total net polarization in a large aperture 
cancels to zero. For inclined or nonaxisymmetric dust 
distributions, that cancellation begins to break down and 
can introduce a systematic bias to underestimate the true 
degree of polarization. In an extreme case, a 100% lin- 
early polarized source would not be measurable with this 
technique at all. But as mentioned above, we restrict our 
efforts to the detection of spatially varying polarization 
from resolved dust. For such systems, the uncertainty 
resulting from this scaling process will be negligible com- 
pared to other factors, such as the biasing effect of the 
PSF halo in Stokes /. 

With that goal in mind — the detection of faint spatially 



resolved polarization from extended circumstellar dust, 
and not the precise measurement of total polarization — 
an IDL data reduction pipeline for IRCAL polarimetry 
data has been developed. This software is available on- 
line^ and automatically handles the reduction, calibra- 
tion, and analysis of IRCAL polarimetry data. Minimal 
user intervention is required to produce flrst-pass science 
reductions, though extensive options exist for overriding 
or fine-tuning the pipeline's behavior. 

6.2. Initial Reduction 

Images are first dark-subtracted, then divided by the 
polarization mode flat field. Sky subtraction is per- 
formed, as described above. Bad pixels are identified 
and removed, using the combination of a mask indicat- 
ing known permanently-bad pixels and the NIRSPEC 
f ixpix iterative bad pixel cleaning routine developed by 
Tom Murphy. At this point, the images are relatively 
clean, and can be corrected for anamorphism using the 
ircal_dewarp routine ( §4.2p . 

Saturated pixels and/or pixels above the detector's lin- 
ear regime are identified by comparison with a threshold 
value of 280,000 electrons pixel"^. The saturated pixel 
lists for the two polarized subimages are merged to en- 
sure that if any pixel is marked as saturated, its corre- 
sponding pixel in the other polarization is also marked 
as saturated. Without this step, polarization artifacts 
would be induced later during image subtraction. 

6.3. Registration 

Accurate registration to a fraction of a pixel is re- 
quired to extract polarized signals very close to a star. 
lApai et al] ()2004l ) found that for AO polarimetry at the 
Very Large Telescope, random shifts of 0.2-0.5 pixel (5-14 
mas) distort the polarized signal of the inner arcsecond, 
and shifts of 1 pixel completely disrupt it. Long exposure 
images, which probe polarization on slightly larger scales, 
are more robust against misalignments. Apai et al. used 
a "two-level Gaussi a n fittin g procedure" to register their 
data, while iPotteil ()2003aD used an iterative shift-and- 
subtract process to empirically find the shift for each 
image that gave the lowest residuals. For IRCAL we use 
Fourier cross-correlation to subpixel precision. Which of 
these techniques is best will vary depending on instru- 
ment and observational conditions. 

The first step in image registration is to determine 
the offset between the two polarized subimages on the 
detector, which varies slightly from night to night due 
to imprecision in the filter wheel rotation. To measure 
this, the full stack of images is summed, and the result- 
ing left- and right-hand subimages are cross-correlated 
to find their displacement. 

The individual images are then registered via cross- 
correlation, fit by Gaussians to subpixel precision. Only 
the left-hand subimages are correlated, with the same 
shifts (plus the beam displacement) assumed correct for 
the right hand field. If necessary, manual intervention is 
possible to tweak the automatic registration. 

An important caveat about IRCAL polarimetry is that 
the images are only Nyquist sampled at Kg. At shorter 
wavelengths, high spatial frequency PSF structure will 

http: //astro .ucla. edu/~mperrin/ircal/ 
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alias to lower spatial frequencies, biasing the correla- 
tions and derived shifts in an unpredictable and time- 
variable manner. In the absence of a priori PSF knowl- 
edge, this imposes a fundamental limit on our abil- 
ity to accurately register these images. Furthermore, 
the long exposures necessary to detect faint circumstel- 
lar dust are frequently saturated on the central star, 
which limits registration accuracy even more. Cross- 
correlation still works to some extent in this case, us- 
ing the diffraction spikes and other outer PSF structure, 
which is why we chose that approach for the IRCAL 
pipeline. The highest accuracy polarimetric image regis- 
tration may require the use of "satellite PS F" techniques 
developed for coronagraphic observatioiis ( | Marois et al.1 
l2006t ISivaramakrishnan fc Oppenheimeill2006f) . 

6.4. Mosaicing 

All the subimages for a given polarization are then mo- 
saiced together, combining different exposure times via a 
weighted average. Subpixel shifting is done by Fourier in- 
terpolation. A separate mosaiced image Mi is produced 
for each of four polarization angles: Mq, M45, M22.5, and 
-^67.5- These mosaics combine polarized subimages from 
both sides of the detector, and images from the redun- 
dant waveplate settings 90, 135, 112.5 and 167.5 are 
summed with their counterparts (i.e. images taken with 
the waveplate at 90° are mosaiced in with those taken at 
0, and so on.) 

The resulting mosaics are then scaled to minimize their 
subtracted residuals Q and U, according to 

Q = Mo - aM45+P 

U = M22.5--fMer.5 + 5 

I = (Mo + aM45 + M22.5 + 7^^67.5) /2. 

The coefficients for minimizing the total power in \Q\ 
and \U\ are determined by IDL's Powell nonlinear opti- 
mization routine. Typical values for a and 7, the scaling 
factors that remove the net polarization of the targets, 
are within a few percent of unity. The (3 and S factors 
correct for any residuals in the sky background levels, 
and are always small. 

This order of operations — mosaicing the images, then 
performing the scaled subtraction — is by no means the 
only way to do things. One can also first difference 
the images and then mosaic the subtracted frames, or 
there are also reduction algorithms that rely upon ra- 
tios rather than differences of the two channels (see 
lTinbergen|[l996[ ). We chose the mosaic-then-subtract 
approach because it offered greater flexibility for deal- 
ing with missing frames from AO system dropouts or 
moments of bad seeing, but the merits of the various 
algorithms should be reevaluated for future instruments. 

The polarized intensity P is computed from Q and U 
in the usual manner, P — \/ Q"^ + U"^- If absolute ac- 
curacy in P were our goal, we would need at this point 
to remove the bias introduced into P by its positive def- 
inite nature, by following on e of the the standard sta- 
tistical debiasing algorithnis (IWardle fc Kronberg|[r974 
ISimmons fc Stewart! [1985 IStewartj|1991[ ). But we skip 
this step for two reasons. First, these corrections are im- 
portant primarily when the signal to noise ratio is low 
{P/ap < 2), and we focus our analysis only on pixels 
with strong detections of polarization {P/<Jp > 3). Much 



more importantly, for AO data the bias in degree of po- 
larization P/I from the uncorrected PSF halo in Stokes 
/ will be far larger than this statistical bias. Recall our 
goal is not to make absolute polarization measurements; 
when attempting to use polarimetry to detect fainter sig- 
nals hidden in the speckled PSF halo, it is simply infea- 
sible to make absolute measurements of P/I and so we 
do not try. 

Along with the mosaics, we generate exposure maps 
and per-pixel empirical uncertainty images for /, Q and 
U, computed from the standard deviation of double- 
differenced pairs of input frames. Computing the 
per-pixel uncertainty based on unsubtracted or single- 
subtracted images would not provide an accurate as- 
sessment of the speckle-noise-reduced characteristics of 
the final Stokes images. During the course of our AO 
polarimetry survey of Herbig Ae/Be stars, these uncer- 
tainty maps proved very valuable in assessing the signif- 
icance of faint signals. 

6.5. Final Steps 

The last steps in reduction are updating the World 
Coordinate System header, and flux calibration in Jy 
arcsec"^ based on the closest-in-time photometric stan- 
dard. The final images are stored as multi-extension 
FITS files containing the Stokes data cubes, exposure 
maps, and uncertainty images. 

These data may be visualized using several IDL plot- 
ting routines, or interactivel y explored u sing a version of 
the atv image display tool ( Barth|[200lD modified to ac- 
cept Stokes data cubes and overplot polarization vectors. 
This code is available from the author's web site^. 

7. ACHIEVED PERFORMANCE 
7.1. Contrast Enhancement 

In practice the IRCAL polarimeter and data reduction 
pipeline remove 98-99% of unpolarized starlight, greatly 
attenuating the background against which scattered light 
must be detected (Figure [3]). The exact level of perfor- 
mance varies with conditions. 

Outside of ^1", we typically reach a noise floor Pnoise 
set by photon and read noise. Inside that radius, instru- 
mental systematics raise the polarized surface brightness 
background above the theoretical noise floor. The precise 
transition between photon-noise-limited and systematic- 
limited depends on factors including source brightness 
and seeing; for some of our targets the transition occurs 
inside 1" (though rarely inside 0.7") and for some it's 
furthe r ou t . This performance is similar to that achieved 
by Po tted (j20033) with the Hokupa'a polarimeter. 

Like speckle noise, the systematic noise floor is not 
Gaussian and does not improve with time as ^/t. Ad- 
ditional integration increases the sensitivity at small 
radii only slightly if at all, similar to the behav- 
ior of se mi-static speckles in conventional AO imag- 
ing (e.g. iRacine et al.l 119991 : lAime fc Soummeii 120041 : 
iHinklev et al.ll2007D . But we stress here that the polari- 
metric noise floor is not merely semi-static speckles rear- 
ing their ugly heads: a perfect polarimeter would com- 
pletely reject any static speckles, no matter how aber- 
rated the / PSF. It is imperfections in the polarimeter 

^ http: //astro .ucla. edu/~mperrin/idl/ 
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Fig. 3. — H band imaging polarimetry of the unpolarized standard HD 18803, demonstrating the cancellation of unpolarized starlight. 
Left: Total and polarized intensity. The display scale for polarized intensity is reduced by 100 compared with total intensity. Right, top: 
Radial profiles showing total and polarized intensity and the minimum theoretical noise floor for P based on photon statistics and read 
noise. As a rough estimate of the light-suppression capabilities of the polarimeter, the thin dashed line shows a copy of the I profile that 
has been multiplicatively scaled by 0.015 to match the P profile. In other words, the polarimeter rejects 98.5% of unpolarized starlight 
for this observation. Right, botto m: A histogram of polarization position vectors relative to the expected centrosymmetric pattern that 
indicates circumstellar dust (sec i|7.2l l. The uniform distribution seen here is the expected null result for an unpolarized star. 

0.04 L J laboratory tests with high-quality polarized calibration 

; + + ; sources. More fundamentally, the various terms in the 

: +^ : error budget all have similar functional forms, Pnoise oc 

■§ 1 /. Flat field errors, registration errors, pupil shear, all 

'■ ^ I introduce noise proportional to total intensity, rendering 

it hard to empirically disentangle the magnitude of each 
individual effect. 

We thus present only a partial resolution to this issue. 
For each target observed during our Herbig Ae/Be sur- 
vey, we estimated the polarized noise floor as a fraction of 
total intensity, by computing the median ratio of the ra- 
dial profiles of P and / between 0.2" and 2.0." In Figure 
H we plot this ratio against the PSF FWHM, a measure 
of wavefront quality. The result is striking: as AO per- 
formance improves (i.e. the FWHM becomes smaller), 
on average a larger fraction of light "leaks through" the 
polarimeter. Better seeing and higher Strehl PSFs per- 
versely result in slightly worse PSF subtractions! 

From this we conclude that the subtracton residuals in- 
crease when the PSF has more power at high spatial fre- 
quencies. Errors in flat fielding would not cause this; such 
errors should be a simple multiplicative factor on the 
PSF, regardless of PSF shape. On the other hand, errors 
from PSF registration, aliasing, and numerical interpola- 
tion will grow as high spatial frequency power increases, 
particularly for an undersampled detector such as IR- 
CAL's. Registration errors of saturated images will also 
contribute to this problem. Long exposures are needed 
to detect faint outer nebulae, and in such data the core 
of the PSF will be saturated and thus unusable for regis- 
tration. In conditions with better seeing, saturation will 
occur earlier, contributing to the observed degradation. 
Thus for IRCAL it appears that image registration and 
aliasing are a more significant limitation than flat field- 
ing. 



Fig. 4. — The fraction of residual starlight not suppressed by 
polarimetry (measured as the ratio of apparent polarized intensity 
over total starlight) plotted versus PSF full-width at half maximum 
(FWHM). The dashed line shows the average suppression ratio in 
each 0.1" wide bin, and the vertical dotted line shows the minimum 
diffraction-limited FWHM. Higher ratios mean more unpolarized 
starlight is leaking through into the final polarized intensity images. 
The suppression of unpolarized starlight is best for observations 
with large FWHM, and conversely worst when the AO correction 
provides the best FWHMs. This surprising result indicates that 
high spatial frequencies in the PSF contribute strongly to the sub- 
traction residuals, suggesting that registration and aliasing rather 
than fiat fielding errors set the greatest limits on performance for 
our polarimeter. 

that allow a fraction of those speckles to bleed from / 
into P and ultimately limit our sensitivity. 

Many factors contribute to this imperfect subtraction 
of images. Is it possible to disentangle which of these — 
flat fielding error, pupil shear, registration error, and so 
on — is dominant? In practice, for IRCAL it has proven 
difficult to quantitatively assess the relative importance 
of the various error terms. In part this is because of 
the instrument's history. It was a mode retrofitted into 
the camera at the telescope, and it never underwent any 
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7.2. Detection of polarized light in reduced images 

In some cases polarized structures are immediately vis- 
ible in reduced images, but more often we face the task 
of detecting faint polarizations at low signal to noise. 
How may we best decide whether a given data set shows 
robust evidence for detected polarization or not? 

Merely looking by eye at polarization images is an in- 
effective way to approach these data. Polarization is fun- 
damentally a vector quantity, and examining only one of 
Q or J7 at a time does not take full advantage of that 
richness. Hunting for a faint signal in P is also not the 
best approach: P incorporates both Stokes Q and U, but 
its positive definite nature biases it so any noise whatso- 
ever results in a nonzero, positive value, which swamps 
faint polarized signals. Numerical fits should be done di- 
rectly on the Q and U images, and only transformed to 
P afterwards if necessary. 

Luckily, light scattered from circumstellar dust has a 
distinct signature we can search for. This characteristic 
may be visualized either as the so-called "butterfly" pat- 
tern in the Q and U images, or as a "centrosymmetric" 
circular pattern of polarization vector s arranged c i rcum - 
ferentially around the central source. lApai et all (|2004f ) 
fit a sinusoid to the butterfly pattern in Q and U images 
to detected scattered light to within 0.1" of TW Hya. 
A closely related technique uses instead the polarization 
position angle 9. The position angles for centrosymmet- 
ric scattering are given by 

dcentro{x,y) = arctau ( - — — ] , 
\x - x^, J 

where {Xi,,y.,) is the central star's location. By com- 
paring the observed position angle 0obs{x,y) with the 
expected position angle centra, we can sensitively dis- 
criminate dust-scattered light from other contaminants 
such as residual speckle noise or registration errors. A 
key advantage is that the position angle Oobs is formed 
from the ratio of Q and U, and hence is less sensitive to 
flat fielding er rors than either of those quantities alone 
(|Potteil |2003ah. 

See Figure [5] for a demonstration of using this method 
to detect a very faint polarized signal from a disk around 
the Herbig Ae star HD 141569. Located at 99 pc, this 
young star has a circumstellar disk believed to be transi- 
tioning from optically thick to optically thin, possibly 
due to the clearing effects of unsee n planets. Based 
on p r ior coronagraphic ob s ervations ([We inberge r et al.l 
119991: IClampin et all 120031 iBoccaletti eTal. .200^ . the 
disk's iJ-band surface brightness is 16-17 mag. arcsec"^ 
between l-2"radii, compared to a stellar magnitude of 
H = 6.8. We observed HD 141569 on 2003 June 17 for 
1600 s total, as part of a large survey of dust around 
Herbig Ae/Be stars (Perrin et al; in preparation). The 
disk's structure is not at all apparent in the / or P im- 
ages. However, if we plot a histogram of the position 
angles for pixels* which have polarized S/N > 3, we find 
they are predominantly close to the expected centrosym- 
metric angles (seen in the lower right panel in Figure [5|) . 
The observed polarization signal is strongest near the 

® In addition to masking out pixels with low S/N, we also mask 
out the diffraction spikes fixed at 45°when computing those cen- 
trosymmetry histograms. 



location of the bright outer ring previously seen by NIC- 
MOS and ACS. While this observation only marginally 
detects this disk and does not reveal its structure clearly, 
to our knowledge this is the first non-coronagraphic de- 
tection of HD 141569's disk in the near-IR. This demon- 
strates that centrosymmetry histograms can be used to 
sensitively detect extremely faint polarized signals which 

might othe rwise be misse d. 

Recently, ISchmid et al.l ()2006[ ) have advocated the use 
of "radial Stokes parameters" as an alternative technique 
that combines aspects of both the butterfiy and cen- 
trosymmetry methods. In this approach the Stokes vec- 
tors are transformed according to 

Qr — +Q cos 2(1) + U sin 2(/) 

Ur = -Q sin 2(1) + U cos 2(j) 

where (j) is the position angle measured from the central 
source. Qr is then a measure of the polarization's cen- 
trosymmetry. The use of these quantities for polarimetric 
analysis deserves further study. 

8. CONCLUSIONS 

We have developed a dual-channel imaging polarimetry 
mode for the Lick Observatory adaptive optics system 
and its IRCAL science camera. Its unique YLF Wol- 
laston prism enables differential polarimetry throughout 
the near infrared, enhancing contrast for the detection of 
faint circumstellar dust. In practice this system can re- 
duce the stellar PSF halo by a factor of 50-100, ultimately 
limited by instrumental systematics such as registration 
and flat-fielding errors. These are the same systemat- 
ics whic h have limited other AO polar imeters such as a t 
Gemini (jPotteHl2003bl ) and the VLT (jApai et al.ll2004f) . 
and we report similar achieved performance to those in- 
struments. 

After more than five years of operations, the IRCAL 
polarimeter has observed targets ranging from young 
stars to planets to planetary nebulae. A few observations 
serve to show the capabilities of the IRCAL polarimeter 
in a variety of scientific contexts. Figure [6] presents the 
first-ever H band imaging polarimetry of Uranus. Fig- 
ure [7] shows high-contrast, multiwavelength polarimetry 
of the circumstellar material around a pre-planetary neb- 
ula. Lastly, Figure[8]presents polarization measurements 
for the synchrotron emission from the inner region of the 
Crab Nebula, demonstrating that the IRCAL polarime- 
ter can also be used for studies of polarizing processes 
besides dust scattering. 

While the IRCAL polarimeter performs well, it suf- 
fers from certain instrumental limitations that could have 
been avoided had polarimetry been considered during the 
initial design of the camera. All too many polarimeters 
are born this way, as post-facto retr ofits to instru ments 
originally intended as pure imagers (jHoughl 120051 ). This 
lesson is particularly relevant to current design efforts 
for "Extreme AO" high contrast imagers: polarimetric 
performance should be considered in trade studies early 
in the design process of both AO systems and science 
cameras. Present AO polarimeters (in every instance 
retrofitted into existing instruments) have ably demon- 
strated their ability to detect faint scattered light in the 
presence of a bright stellar PSF — the potential for dis- 
covery using a truly optimized AO polarimeter is tremen- 
dous. In IRCAL's case, one of the major difficulties is ac- 
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Fig. 5. — H band polarimetry of a very faint signal: the circumstellar disk of the Herbig Ae star HD 141569. See the caption of 
Figure |3] for a description of the various quantities plotted. The dotted ellipse traces the outer bright ring in the disk, as seen by ACS 
coronagraphy ( Clampin et al. 2003). In these IRCAL data, little to no sign of circumstellar dust scattering is apparent in the images or 
the radial profiles. However, the histogram of position vectors has a definite peak near zero, indicating that the majority of pixels above 
the polarized S/N threshhold have centrosymmetric polarization angles, as expected for dust scattering. While these data do not provide 
much scientific insight, this is a statistically significant detection of polarized light, demonstrating the ability of differential polarimetry to 
sensitively detect very faint signals for subsequent followup. 

registration of images with saturated (or coronagraphi- 
cally occulted) PSF cores. 

An even more promising technique would be to split 
the polarizations only after pixellating the image plane, 
entirely eliminating the subpixel registration problem 
and partially mitigating the effects of uncertain flat 
fields. Several recent designs for high contrast imaging 
spectrographs use focal plane lenslet arrays to chop the 
PSF into pixels prior to wavelength d ispersion, to reduce 
non-c ommon-path optical errors (e.g. iMarois et ai]|2004 
|2005( ). Experiments to validate these designs a re cur 
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rently underway both in the laboratory ( Lafreniere et al 
2006) and on the sky with OSIRIS at Keck (jLarkin et al 



20061). Recent results are promising, but differential re- 
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Fig. 6. — Uranus, as seen in polarized H band light on 2004 
July 05. The upper frame shows total intensity while the lower 
frame shows polarized intensity with overplotted vectors giving 
the percent polarization. Uranus is very limb-brightened in po- 
larized light, with a peak polarization of about 1.5% at the edge of 
its disk. Unlike circumstellar disks, the polarization vectors here 
point radially inwards rather than circumferentially, an expected 
effect of Rayleigh scattering in planetary atmospheres. The north- 
ern hemisphere, at right, is tilted away from us and noticably less 
polarized than the nearer southern hemisphere. To our knowledge, 
this is the first resolved imaging polarimetry of Uranus in the near- 
infrared, and only the second imaging pola rimetry of the planet at 
any wavelength, after ISchmid et al.l 1120061) . 

curately registering undersampled or core-saturated im- 
ages. Future polarimeters should avoid this by having 
a well-sampled detector and using "satellite PSFs" , es- 
sentially intentional ghost images, to allow the precise 



fraction (both atmospheric and in AO system dichroics) 
poses complications. However, for polarimetric applica- 
tions, differential refraction is much less of an issue, and 
the ability to modulate the polarization provides robust- 
ness against non-common-path errors. Lenslet-based dif- 
ferential polarimeters deserve careful study as part of the 
high contrast astronomer's future toolkit. As an added 
bonus, such an instrument would be immune to lateral 
chromatism, eliminating the need for exotic materials 
such as YLF. 

We end with a brief list of lessons learned from IR- 
CAL, applicable to the design of future high contrast 
polarimeters. Starlight suppression of about 2 orders of 
magnitude is attainable with IRCAL and all its flaws. 
A goal of 3 orders of magnitude suppression seems not 
unrealistic for future instruments. 



Low repeatability of the aperture and filter wheel 
mechanisms in IRCAL complicates data taking and 
reduction, but is not insurmountable. Any reason- 
able future design will improve on IRCAL in this 
regard. 
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Fig. 7. — Multiwavelength imaging polarimetry of the pre-planetary nebula Hen 3-1475. This data shows the power of polarimetry to 
reveal faint circumstellar material even in the presence of complicated and variable AO PSFs. In total intensity (top row) the circumstellar 
nebula is not easily visible, but in polarized intensity (bottom row) the symmetric bipolar nebula is clearly d etected in all three wavelengths. 
The morphology seen here compares well with recent HST NICMOS observations presented bv lUeta et al.l (.2007.) . 




Fig. 8. — Laser guide star AO imaging polarimetry of the inner portion of the Crab Nebula, Ml, in H band. AO polarimetry can be 
applied to observe sources of polarization besides dust scattering. Here we observe near-IR synchrotron emission from the Crab Nebula. 
The arrow indicates the Crab's central pulsar. Several well known wisps and the knot immediately southeast of the pulsar are visible. The 
observed morpholo gy and polari zation angles are generally consistent with the higher angular resolution observations obtained recently 
with HST ACS bv IHesteH 1)20071 ) . Measurement of the degree and orientation of polarization can provide insight into the magnetic field 
geometry and relatavistic outflow properties. The large field of view shown here was obtained by mosaicing together dithered images 
comprising 5700 s total exposure time, taken on 2005 Nov 21. 



2. Repeatability of mechanisms is particularly impor- 
tant as there are stringent requirements on flat 
fielding, and IRCAL's mechanism nonrepeatabil- 
ity prevents the development of high signal-to-noise 
flat libraries. 

3. The WoUaston prism should ideally be located after 
the cold pupil, to prevent pupil shear between the 
two polarized beams. 

4. To minimize instrumental polarization, the wave- 
plate should be located as far upstream as possi- 
ble. However, optimized scaled subtraction offers 
an effective way to empirically remove instrumen- 



tal (and interstellar) polarization for applications 
where absolute polarimetry is not required. 

5. Registration errors are a key component of the er- 
ror budget, particularly for undersampled data or 
long exposures saturated on the central star. The 
use of artifical "satellite PSFs" as astrometric ref- 
erences may help solve this problem. 

6. Lenslet arrays hold great promise for future high 
contrast polarimeters, largely eliminating the im- 
age registration problem. Post lenslet array, the 
requirements on lateral chromatism are tremen- 
dously relaxed, potentially greatly easing require- 
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ments on the beamsplitter. 

7. For detecting very faint polarized signals, cen- 
trosymmetry histograms provide a useful detection 
metric, particularly when coupled with numerical 
uncertainty maps computed as part of the data re- 
duction process. 
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APPENDIX 

AN UPDATED LIST OF MATERIALS FOR INFRARED WOLLASTON PRISMS 

A la rge number of birefringent optical materials can be used to make astronomical WoUaston prisms. lOliva et al.l 
(|1997( ) provides a useful list of some of the best candidates. Since that paper, updated measurements of the optical 
properties of many of these materials have become available, while other candidate materials have come to light. Table 
[1] lists potential materials for near infrared WoUaston prisms, and provides updated values for the birefringence An 
and birefringent dispersion parameter Vau, while Figure [9] shows both achievable beam separations and birefringent 
optical properties for there materials. 




Fig. 9. — Left: Beam separation angle as a function of WoUaston prism angle, for birefringent materials from Table [Tl More birefringent 
materials are higher in this plot, and require a smaller prism angle to acheive a given beam separation. The horizontal line marks the 
desired separation for IRCAL. Right: Birefringence versus birefringent dispersion for these materials. This is the birefringent version 
of the standard Abbe plot for optical materials. On the horizontal axis, the birefringe nt d ispersion parameter Vau measures how much 
birefringence varies with wavelength, with higher numbers indicating less variation (see il3.2| l. Imaging polarimeters benefit from materials 
whose birefringence is more constant with wavelength (high VAn), such as YLF and AgGaS2- For each material, we show the change in 
these birefringent parameters from 1-2.5 fj,m. Triangles, circles, and squares denote the central wavelengths of J, H, and Kg respectively. 
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TABLE 1 

Materials for Infrared Wollastons 



Material n An VAn,J VAn,H VAm K References 



GaSe 


2.74 


0.3195 


157 


135 


484 


1 


Ag3AS3 


Z. ( 1 




A 1 


( z 


lUo 


z 


YVO4 


2.15 


0.2053 


97 


172 


278 


3 


CaCOs 


1.63 


0.1550 


50 


33 


21 


4 


Ag3SBs3 


2.83 


0.1424 


37 


66 


102 


2 


Li03 


1.85 


0.1367 


107 


109 


86 


2 


BBO 


1.64 


0.1158 


1367 


175 


87 


2 


LiNbOa 


2.21 


0.0728 


48 


45 


33 


2 


AgGaS2 


2.42 


0.0535 


227 


4152 


11210 


5, 6 


BABF 


1.61 


0.0397 


19 


14 


9 


7 


ADP 


1.48 


0.0261 


7 






8 


YLF 


1.44 


-0.0220 


268 


1072 


169 


9 


KDP 


1.48 


0.0208 


5 






10 


MgF2 


1.37 


-0.0113 


135 


89 


56 


11 


AI2O3 


1.74 


0.0078 


235 


213 


122 


12 


KTiOAs04 


1.79 


-0.0055 


56 


26 


17 


1 



Ref erence s. — 1: [A llaklivcrd iey et al.l II2005I') 2: Weberl 
||2003') 3: ILomheim fc DeS hazer (197^) 4: IGhoshI II1999I1 5 : 
IWiller et aLl J2OOII) 6: Takaoka & Kato ( 1999) 7: H u et al.l (1200^ ') 
8; CASIX data sheet, http: / /www.u-oplaz.com/crystals/c rystals07.htm| 
9: [Barnes & Gettcmy ( 198UI) 10: Kedoptronics Uafa sheet, 
|http: / /www. optical-componcn ts.com /KDP-crystal. html 11: ITropfl 
II19!J5I) 12: IKaplan Hz I'homasI ( 12003) 

Note. — Materials are ordered by decreasing birefringence. The index 
of refraction, n, and birefringence, An, are stated for 1.65 iim. The 
polarization dispersion parameters Vs„ were calculated assuming filter 
bandpasses of 20%, centered on 1.25, 1.65, and 2.1 ^m. Dashes indicate 
wavelength ranges where accurate index of refraction information was 
not available. More detailed tables of optical properties and software 
implementing Sellmeier models for all these materials is available from 
the authors. 
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